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The  backscatterlng  of  light  from  drop lots  of  motor  in  clouds  is  known 

to  bo  enhanced  duo  to  a  weak  focusing  In  the  backward  axial  direction.^  This 

axial  focusing  gives  rise  to  the  optical  effect  known  aa  the  "glory"  (Kef.  1-6 

and  references  cltsd  therein).  Axial  focusing  Is  also  known  to  be  present  in 

the  clssslcal  and  quantum  mechanical  scattering  of  particles  where  the  enhance- 

ment  of  the  differential  cross  section  Is  known  as  the^  "glory  effect."  Mani- 

8  9 

festatlons  of  this  type  of  focusing  In  acoustics.,  though  previously  noted,  * 
are  relatively  unexplored.  In  the  present  paper  we  derive  a  physical-optics 
approximation  which  describes  '.he  effects  of  diffraction  on  the  axial  focusing 
of  glory  rays  In  fluid  sphares.  He  verify  our  approximation  by'  coopering  It 
with  numerical  computations  of  the  partial-wave  sum  for  the  exact  scattering 
from  lnvladd  fluid  spheres  (formulated  e.g.  in  Kef.  10).  It  Is  also  shown 
that  the  backscattarlng  can  be  further  enhanced  by  choosing  the  sound  velocity 
ratio  such  that  glory  and  rainbow  rays  coincide. 

The  paths  of  rays  through  a  sphere  are  determined  by  the  acoustic 
refractive  index  H  •  c  /c,  where  c  and  c„  are  the  sound  velocities  of 

O  1  O  X 

the  outer  end  Inner  fluids,  respectively.  The  Incident  wave  Is  taken  to  be  a 
plane  wave;  It  is  unmodulated  and  has  a  wavelength  X  In  the  outer  fluid.  A 
physical-optics  approximation  for  the  diffraction  limited  axial  focusing  Is 
derived  here;  the  derivation  assumes  that  X  «  the  sphere's  radius  a. 
Comparisons  with  the  exact  partial-wave  sum,  made  for  several  M  <  1  with 

ka  •  100  and  1000  (where  k  ■  2rr/y  i«  the  wavenumber  of  the  Incident  wave) , 
demonstrate  the  legitimacy  of  the  method.  Kays  which  reflect  from  the  front  end 
rear  poles  of  the  sphere  (the  "axial  rays")  experience  no  focusing;  consequently 
their  amplitudes  are  smaller  than  those  of  the  focused  glory  rays  when  ka  Is  large. 
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A  second  paper ^  extends  the  present  results  to  csrtsla  esses  of  glory 
scattering  of  ultrasonic  pulses  by  an  elastic  sphere  in  water;  it  also 
describes  direct  observations  of  diffraction  limited  backward  axial  focusing. 

When  rays  inside  an  elastic  sphere  are  not  mixed  la  their  type  (l.e. ,  they  are 
all  shear  or  all  longitudinal  rays),  ths  paths  are  the  sane  as  those  for  a 
fluid  sphere  with  c^  taken  to  be  either  the  shear  or  the  longitudinal  wave 
velocities.  Consequently,  the  paths  in  fluid  spheres  with  M  <  1  most  closely 
resemble  those  for  most  elastic  spheres.  The  emphasis  of  the  present  paper  is  on 
fluid  spheres  with  H  <  1;  however,  in  Sec.  7  we  model  scattering  from  fluid  spheres 
with  M  >  1  and  present  computations  of  the  exact  scattering. 

The  present  description  of  the  acoustic  glory  may  also  be 
extended  to  the  case  of  a  spherical  elastic  shell  filled  with  a  liquid.  Targets 
of  this  type  are  asserted  to  be  useful  as  navigational  aids  and  for  the  calibration 
of  sonar  devices.  As  described  in  Sec.  7,  backs  cat  ter ing  from  a  fluid  sphere 
should  be  quite  large  when  M  *  1.180.  This  anhancenent  of  the  scattering, 
which  is  due  to  a  coincidence  of  rainbow  and  glory  rays,  should  also  be 
applicable  to  the  design  of  liquid-filled  shells  with  unusually  large  target 
strengths.  Some  applications  to  shells  are  noted  in  Sec.  71. 

Our  derivation  of  the  scattering  amplitude  for  backward  and  near  back¬ 
ward  directions  will  parallel  our  previous  treatmsnts  of  the  backseat terlng  of 

4*6  5 

light  from  spherical  bubbles  in  liquids  and  other  dielectric  spheres.  To 

facilitate  a  comparison  with  the  exact  scattering,  we  derive,  for  the  first 

time,  the  dependence  of  the  phase  of  the  glory  scattering  on  the  scattering 

eagle. 

1.  Bay  Acoustics  sad  Axial  Tosuslag 

In  this  section  we  review  the  elementary  ray  acoustics  of  a  sphere  aad 
the  geometric  effect  which  gives  rise  to  axial  focusing.  Kay  paths  are 
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described  by  Snell's  law: 


•ln6  -  M  slnv  (1) 

where  6  la  the  angle  of  incidence  et  Che  sphere's  surface  and  v  la  the 
angle  of  refraction.  The  auaber  of  chords  Inside  the  sphere  for  e  given  ray 
will  be  denoted  by  n.  Several  rays  which  are  reflected  and  refracted  la  the 
near  backward  direction-  ere  shown  In  Fig.  1.  The  deviation  Y  In  the  direction 
of  a  ray  fron  the  backward  axis  (the  CC'  axis)  is  given  by 

Y  -  e  -  B  ,  (2) 

6  •  2nv  -  0  +  (2  -  n)ir,  n  >  1  (3) 

where  6  Is  the  Internal  angle  (relative  to  the  CC'  axis)  of  the  point  et 
which  a  ray  leaves  the  sphere.  (Equation  (2)  nay  be  derived  by  noting  that 
the  direction  shifts  due  to  the  initial  and  final  refractions  have  the  sane 
aagnlcude.  See  Fig.  2  .]  The  roots  of  y  •  0  having  0  <  8  <  0^  will  be 
denoted  by  0Q  and  the  associated  rays  will  be  referred  to  as  glory  revs.  Here 
0___  -  ir/2  whan  M  >  1  and  0__  »  0  whan  M  <  1  where  0  ■  arcsln(M)  is 

MX  MX  C  C 

the  critical  angle  of  incidence.  Exact  backscatterlng  also  occurs  for  rays 
with  8*0  and  Y  “  (n-2)ir,  n  •  0,2,4  .  .  .j  these  rays  will  be  referred  to 
as  axial  revs.  Certain  axial  and  glory  rays  are  shown  as  the  solid  lines  In 
Fig.  l  . 

The  following  la  a  suanary  of  conditions  on  N  and  n  for  glory 

12 

rays  to  exist.  It  is  well  known  that  n  •  2  glory  rays  exist  for 
/l  <  M  <  2.  The  eaph asls  of  the  present  paper  will  be  e  dees  of  glory  rays 
which  exist  for*  0  <  Ki  with  a  >  2.  The  upper  bounds  are  >  1 
with  ■*  1  as  a  •  and  We  have  previously  demonstrated 

that:* 
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(31/2  6  -  9) 1/2  a  1.179960 

(4) 

(4/3)(2/3)1/2  «  1.088662 

(5) 

Moat  of  our  discussion  vlll  be  concerned  with  spheres  with  M  <  1. 

The  general  class  of  glory  rays  in  such  spheres  Is  described  by  0  »  8' 
where  6*  »  0  +  n'2ir  where  2n'  +  1  Is  the  number  of  times  rays  cross  the 
symmetry  axis.  It  Is  necessary  for  n  >  2n'  +  2,  n'  ■  0,1,2,  .  .  .  and  conse¬ 
quently  glory  rays  with  n*  >  0  have  4  or  aore  Internal  reflections.  Ve  have 
extended  the  physical-optics  approximation  described  in  Sec.  II  and  III  to 
Include  rays  with  n*  >  0  and  find  that  their  contribution  to  the  total 
glory  scattering  from  fluid  spheres  is  significantly  smaller  than  that  of  rays 
with  n’  -  0.  Consequently,  for  reasons  of  brevity  and  simplified  notation, 
we  will  only  describe  the  glory  scattering  due  to  rays  with  n*  •  0  which  is 
the  class  of  rays  illustrated  in  Fig.  1.  Irrespective  of  the  value  of  n',  as 

n  -*■  «,  0  -*•  0  . 

Q  C 

Consider  rays  which  lie  close  to  a  glory  ray;  e.g. ,  the  dashed  lines 
near  the  n  -  3  glory  ray  In  Fig.  1.  The  incident  wavefront  bounded  by  these 
rays  is  an  annular  ring  which  corresponds  to  wavelet  de  in  Fig.  1.  The 
corresponding  wavefront  Is  toroidal  when  It  leaves  the  sphere  and  It  corres¬ 
ponds  to  curved  wavelet  d'e’ .  This  wavelet  is  toroidal  because  the  figure  nay 
be  rotated  around  the  CC*  axis.  The  nth  toroidal  wavefront  appears  to 
originate  at  a  virtual  rlagllka  source  at  F  .  Each  source  forms  a  virtual 

Q 

focal  circle.  Bays  la  the  portion  of  the  outgoing  wavelet  F'd'  cross  the 
backward  axis  when  they  are  extended. 

When  the  scattering  amplitude  from  a  penetrable  sphere  is  computed  via 


ray  optics,  the  amplitude  diverges  as  the  observation  point  approaches  the 
backward  axis  [see  e.g.,  Bef.  1  (Sec.  12.21),  Bef.  9  (Sec.  IC3),  or  Bef.  13 
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(S«c.  VIII)].  This  divergence  is  the  manifestation  of  the  exiel  focusing  sad 
it  also  effects  the  scattering  along  the  forward  axis .  ®  ®  This  divergence  Is 
also  present  in  the  classical  description  of  scattering  froa  a  central  potential 
where  the  well  known  differential  cross  section  ls^ 

f  -  lib  i  bii'V1»i  <« 

where  the  sua  Is  over  the  different  particle  trajectories  which  are  scattered 

by  an  angle  $  and  is  the  iapact  paraaeter  of  the  1th  class  of  trajectory. 

Axial  focusing  is  predicted  for  those  rays  with  bjJdbj/d^l  i4  0  as  $  ♦  ir 

(backseat taring)  or  $-*•()  (forward  scattering). 

The  cause  of  the  focusing  can  be  seen  by  rotating  Fig.  1  about  the 

CC’  axis.  Follow  adjacent  rays  having  Infinltesiaally  different  aslauthal 

angles  but  having  the  saae  angle  of  Incidence  0.  When  this  0  Is  slightly 

less  than  the  0  of  some  "glory  ray,"  the  adjacent  rays  cross  the  axis  at  a 
n 

common  point  after  they  leave  the  sphere.  (This  can  be  seen  by  extending  the 

line  F^d'.)  This  crossing  of  aslnuthallv  adjacent  rays  gives  rise  to  the 

geometrically  predicted  divergence  of  the  energy  density  on  the  axis. 

The  method  for  correcting  for  this  divergence  was  suggested  by 

Van  de  Hulat.*’1^  It  makes  use  of  a  physical-optics  approxin**"*""^^  which 

Involves  (a)  the  computation  of  amplitudes  in  an  exit  plane  near  the  sphere 

via  ray  optics,  and  (b)  allows  these  waves  to  diffract  to  the  observation 

point.  This  method  will  be  Illustrated  for  glory  rays  which  give  rise  to 

toroidal  wavefronts  such  as  those  shown  In  Fig.  1.  It  will  be  shown  in  Sec.  V, 

however,  that  as  M  ■*>  M' ,  the  nth  wavefront  Is  no  longer  toroidal. 

n 

There  are  rays  and  types  of  scatterers  for  which  there  is  no  geomet¬ 
rically  predicted  focusing.  There  Is  no  focusing  of  backscattered  axial  rays 
except  for  certain  M  >.  2  (e.g. ,  n  ■  2  and  K  •  2,  see  Sef.  1A).  The  n  •  0 
ray  reflects  without  entering  the  sphere.  To  sn  external  observer.  Its 
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reflection  appear*  to  corns  f roa  a  point-like  source  at  Ag  (see  Vlg.  1) .  It 

Is  also  evident  that  when  a  plane  wave  la  incident  on  a  cylinder  (with  M  <  1) 

and  propagating  in  a  direction  perpendicular  to  the  cylinder's  axis,  there 

8  9 

will  be  no  focusing  of  the  backs cattering.  '  Figure  1  Is  applicable  but 
with  rotation  about  the  CC'  axis  no  longer  allowed.  Backscattered  rays 
appear  to  originate  from  virtual  line  sources  located,  e.g.,  at  Aq  and  the 
Fq.  It  Is  evident  fron  asymptotic  formulations  of  scalar  wave  scatter'  by 
cylinders3"** 3-6  that  tte  geometric  scattering  from  a  fluid  cylinder  dost 
diverge  as  y  •*  0. 

II.  Amplitude  and  Phase  In  the  Exit  Plane 

In  this  section  we  use  ray  optics  to  describe  the  amplitude  and  phase 
of  the  glory  waves  In  the  exit  plane.  It  la  convenient  for  this  plane  to  be 
the  one  which  touchee  C*  with  Its  normal  parallel  to  the  propagation  direc¬ 
tion  of  the  Incident  wave.  Its  projection  onto  Fig.  1  Is  the  dashed  vertical 
line.  After  the  incident  ray  crosses  the  dashed  vertical  line,  the  propagation  . 
phase  delay  for  reaching  the  exit  plane  is: 

n  •  k[a(l  -  cos0)  4-  2anM  cosv  +  w]  (7) 

where  (Fig.  2)  w  Is  the  distance  traveled  by  the  ray  from  the  exit  point  on 
the  sphere  to  the  exit  plane: 

w  -  a(l  -  cos6)  sec(8  -  8)  .  (8) 

The  ray  crosses  the  exit  plane  at  a  radius  s  from  C'  where  from  Fig.  2 

s  »  a[slnB  -  (1  -  cosB)  tan (6  -  8))  (9) 

The  radius  of  the  nth  focal  circle  Is  b  ■  s(sln8  ).  As  Is  evident 

u  XL 

from  Fig.  1  and  by  direct  computation,  dn/ds  •  0  when  s  •  b  .  The  radius 

n 


a  of  the  toroidal  wavafront  at  tha  axlt  plana  la  obtainad  by  computing  tha 
n 

wavafront' a  curvatura: 

a  -  k^n/de2)"1  ,  a  -  b  ,  (10a) 

n  n 

-  a[l  +  |(nt  -  l)"1  cos0J  ,  (10b) 

z  n 

whara  t  ■  tanv  /tan6  and  v  la  given  by  (1)  evaluated  at  9  ,  tha  glory 
q  q  q  n 

ray  condition.  Tha  proof  of  Eq.  (10b)  la  outlined  in  Appendix  A.  The  apraadlng 
of  tha  wavelet  da  at  the  exit  plana  la  characterized  by: 


where  (d'a')  and  (da)  are  the  arc  and  linear  lengths  of  the  outgoing  and 
incoming  wavelets,  respectively,  and  s  is  the  distance  from  C'  to  the 
indicated  points  of  contact  with  the  dashed  plane.  From  symmetry  arguments, 
the  rightmost  side  of  (11),  when  generalized  to  arbitrary  n,  becomes: 
q_  •  |lim[b  -  a(0)]/(b  -  a  sin0) I ,  as  0  ■*  0  where  a  is  an  Implicit 
function  of  the  ray's  original  angle  of  incidence  via  (1),  (3),  and  (9). 
Application  of  L'Hoapltal'a  rule  to  this  limit  (see  Appendix  A)  gives: 

-  |1  +  2(nx  -  l)aec0n|  -  |ac/(aQ-a)|.  (12) 

With  M  <  1,  a  >  a  for  all  finite  n  (and  n')  ao  that  the  abaolute  value 

U 

eigne  in  (12)  are  not  needed;  for  rangea  of  M  >  1,  however,  we  find  aome 

(a  /(a  -a)]  are  negative. 

II  s 

Let  Pj  exp(-lut)  denote  the  incident  preaaure  in  the  daahed  plane 
through  C'  in  Fig.  1  where  u  and  p^  are  the  frequency  and  amplitude  of 
the  wave.  Say  optica  givea  the  following  amplitude  in  the  exit  plane  for  the 


lim 
(de)  - 


.(-d'-O  . 

0  (de) 


lim 
(de)  - 


s(e)  -  s(d) 
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nth  toroidal  wave: 

,;<•>  .  ♦%♦»(.-  bp)2/^]  (u) 

Si 

where  the  toroid  haa  been  approximated  by  a  quadratic  surface.  This  approxima¬ 
tion  introduces  a  negligible  phase  error  provided  k(s  -  b  )*  «  la^l.  The 

n  q  1 

phase  factor  nQ  Is  the  propagation  phase  delay  of  the  glory  ray: 


n  ■  2ka(l  -  cos 6  +  nMcosv  )  (14) 

n  no 

and  Is  the  phase  shift  due  to  the  crossing  of  focal  lines. ^  The  q^2^2 
factor  accounts  for  the  change  In  the  area  of  the  wavefront. 

The  factor  B(s)  accounts  for  the  reduction  in  amplitude  due  to  the 
partial  transmission  or  reflection  of  the  wave  at  each  interface.  We  approxi¬ 
mate  this  by  repeated  use  of  the  Internal  reflection  coefficient  for  pressure  at 
*  Plmom  surface^  ^  as  a  function  of  the  external  angle  8 


8(6) 


(M2  -  sln28)1/2  -  TcosQ 
(M2  -  sin20)1/2  +  Tcos6 


(15) 


where  T  “  p./p  and  p.  and  p_  are  the  densities  of  the  Inner  and  outer 

X  O  X  o 

fluids,  respectively.  In  (15),  8  Is  chosen  to  be  angle  of  incidence  of  the 
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ray  which  crosses  the  exit  plane  at  s.  Symmetry  relations *  between  the 
transmission  and  reflection  coefficients  yield  the  following  combined 
coefficient 

B(s)  -  R(0)Bf’1  [1  +  R(6) ] [1  -  R(8) 1  .  (16) 

The  phase  shift  accounts  for  the  phase  advance  of  ir/2  associated 

with  each  crosalng  of  a  focal  curve  prior  to  reaching  the  exit  plane.  (This 

1  19 

shift  occurs  due  to  vanishing  of  the  wavelet's  area  at  each  focal  curve.  ’  ) 


9 


From  Fig.  1  it  is  evident  that  there  ere  two  types  of  curves.  One  type,  et 

points  in  Fig.  1,  is  due  to  the  intersection  of  initially  adjacent  rays 

which  lie  in  the  sene  meridional  plane.  There  are  n  -  1  focal  curvea  of 

this  type  when  M  <  1.  The  second  type  is  due  to  axial  focusing  of  rays 

within  the  sphere.  There  is  one  focus  of  this  type  each  tine  the  internal  ray 

Intersects  the  axis;  this  occurs  once  when  n'  •  0.  This  is  the  point 

for  the  n  ■  3  ray  in  Fig.  1.  The  total  shift  becomes  u  •  -  mr/2.  (This 

n 

result  differs  from  that  given  in  Bef.  1,  Sec.  12.22  as  Van  de  Hulst  was  con¬ 
cerned  with  the  phase  shifts  at  a  distant  observer.) 

Evaluation  of  the  constants  b  ,  a  ,  q_,  and  n  requires  that  (1) 

n  n  a  n 

through  (3)  be  solved  for  8  with  y  ■  0.  Vlth  n  •  3  or  4  the  system  of 

n 

equations  reduces  to  a  cubic  equation  which  leads  to  the  following  result3 


•in2  8  -  h’M^l  -  h  cos  [(I*  +  ir)/3]>  (17) 

n  u  n  n 

r3  -  arccos  [  (1  -  K^h”3]  (18a) 

rA  “  arccos  [(27/16)!*2  -  1]  (18b) 

where  0  i  ^  i  i,  h’  ■  1/2,  h3  -  [1  -  3“1K2I**,  h'  -  2/3,  and  h4  -  -  1/2. 

When  n  >  4,  we  solve  the  system  Iteratively  for  6  by  choosing  8  ,  as  the 

n  n-i 

first  estimate  of  0  . 

n 

111.  Diffraction  and  the  Far-Fleld  Scattering  Amplitude 
A.  Stationary-Phase  Approximation  for  the  Backscattered  Amplitude 


The  amplitude  p  at  the  observation  point  Q  can  be  expressed  in 
n 

terms  of  e  diffraction  integral  of  P^  1b  the  exit  plane  (the  x'y'  plane 

in  Fig.  3).  In  Fig.  3  the  s  axis  is  the  extension  of  the  CC'  axis  and  the 


10 


backscatterina  angle  y  Is  measured  with  respect  to  C*  and  Cba  s  axis. 

7t^  y '  axis  nay  b«  chosan  to  be  Che  daahad  vertical  axis  in  Fig.  1.  (Unlike 

4 

the  corresponding  optics  problem  in  which  polarisation  breaks  the  syanetry, 
here  the  orientation  of  the  y'  axis  is  arbitrary.)  The  Fraunhofer  approxl- 

io  2Q 

nation  *  for  r’  gives  r'  *  r  -  [(xx1  ♦  yy’)/r]  and  it  gives  the 
following  approximation  of  the  diffraction  integral  for  p 

n 

p„  -  lirUAr)”1  eikr  F,  r  »  ks2  ,  (19) 

n  wit 


F(x,y) 


-  £  fl  p;(«1,y).-Ut(,°‘’  + 


(20) 


*•  V2  <  .2 


where  aM<>  the  radius  of  the  domain  of  integration,  is  sufficiently  large 
that  the  contribution  to  the  diffracted  amplitude  from  points  outside  of 
is  negligible.  In  addition  to  making  use  of  the  Fraunhofer  approximation, 
(19)  neglects  any  corrections  to  the  amplitude  due  to  the  obliquity  between 
Q  and  the  wavefront  in  the  exit  plane.  He  anticipate  the  result  below  that 
the  phase  of  the  integrand  ia  stationary  in  the  region  of  applicability  of  the 
approximations  leading  to  (13).  Conaequently  we  extend  *MZ  in  (20)  to  “> 
and  make  use  of  the  circular  symmetry  of  p'  by  writing  F  in  terms  of 
radial  and  azimuthal  integrals2®’22- 


F(Y) 


C 


•  H(s,y)de 


(21) 


i  r2* 

W(e,Y)  -  I  exp(-iker  cos (i|»  -  £)]d*  -  JQ(u)  (22) 

2  2  7  1/2 

where  T  -  [(x  ♦  y  )/r  ]  -  slny,  u  “  ksl\  and  ip  la  the  aslmuthal  angle 

of  (x' ,y').  Due  to  aymMtry ,  F  and  H  do  not  depend  on  the  aslmuthal  angle 
$  of  Q. 
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Sine*  W  Is  given  by  chs  xero-order  Bssssl  function  JQ(u),  F  Is  Chs 

Hsaksl  transform'’ ‘  of  p'.  Ths  phass  of  p'  is  stationary  whan  s  »  b 

no  n 

and  ww  usa  ths  atatlonary  phass  approximation  (SPA)  of  F.  Ths  amplitude  of 
p*  is  proportional  to  B  which  Is  a  slowly  varying  function  of  s  near  b 

a  u 

sxcspt  when  0  is  doss  to  ths  critical  angle  0  .  As  n  •,  recall  that 
0Q-^  ©c»  &  -*■  ±1,  and  B  0.  Consequently,  the  strongest  pQ  are  those  with 
snail  n  where  B  is  sufficiently  slowly  varying  that  it  nay  be  removed  from 
chs  integrand.  These  approxina clone  give 

F  *  Pib<V‘£1/2  ♦  Wn)]  I  (23) 

I  -  j^s  JQ(u)  sxp[ik(s  -  bn)2/2an]ds  (24) 

-1/2 

where  q^  has  also  been  reaoved  fron  the  integrand  since  ths  spreading  of 

the  wavefront  has  been  approximated  by  (12)  which  is  the  value  at  o  ■  b  .  The 

n 

SPA  of  I  is  (see  e.g.  Baf.  22,  Sec.  4.2c): 

1  *  bn(kbn)~1/2  {(2*iaalbn>1/2  W*1"74  +  *  0(kV’1}  (25) 

where  u  •  kb  r  and  the  order  of  the  correction  tern  was  obtained  by  noting 

Q  Q 

that  the  lowest  order  corrections  fron  the  endpoints  of  (24)  vanish.  The  SPA 
requires,  however,  that  J_(u)  be  slowly  varying  near  u  •  u  and  consequently 

V  U 

2 

that  y  be  snail.  The  condition  on  (19)  nay  be  written  r  »  kbQ  since  F 
is  dominated  by  contributions  to  ths  integral  with  s  3  b  . 

B 

B.  Deterainatlon  of  the  Phase  fron  Properties  of  the  Angular  Spectrun 

When  Q  is  not  dose  to  the  backward  axis,  (25)  does  not  yield  the 
correct  phase  for  p_  due  to  the  requirsnent  that  u  be  snail.  The  phase 

u  » 

correction  is  deternieed  hers  by  shifting  the  plane  of  integration  to  that  of 
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the  focal  circle.  (The  point  at  which  tha  axis  crosaaa  this  plana  will  ba 

danotad  by  C  ;  C,  ia  shown  in  Fig.  1.)  Sinca  wa  lat  a  -*•  •  la  (20), 
n  j  wax 

F  is  a  Fourlar  transform  which  expresses  tha  planavava  angular  spaetrun  of 
p'.  Tha  spaetrun  in  tha  salt  plana  is  ralatad  to  tha  spaetrun  F  ,  determined 

U  Q  " 

in  tha  shifted  plana,  of  tha  virtual  source  appropriate  for  tha  shifted  plana. 

21 

The  relation,  which  is  wall  known  la  optica,  is  given  by 


F  -  exp (lka  cosy)F_ 
ii  n 


(26) 


where  a  is  the  distance  between  tha  two  planes  and  ket  cosy  is  the  phase 
q  n 

shift  of  a  planavava,  tilted  at  an  angle  y  with  respect  to  tha  axis,  as  it 
travels  from  the  plana  at  C  to  tha  one  at  C' . 

II 

In  tha  shifted  plana,  tha  virtual  source  is  a  circle  whose  spatial 

properties  are  those  of  tha  radial  6  function:  6(s  -  b  ).  Writing  F  in 

n  n 

tha  form  of  a  Hankal  transform  gives: 

Fo  ■  V">  «<•  -  V*  ■  'AW  <27> 

where  Dq  is  a  complex  constant  which  has  bean  determined  by  requiring  DnbnJg(un) 
•xp(ikotQ)  •  F  via  (23).  This  procedure,  combined  with  (26)  and  (19),  gives: 
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where  (31b)  makes  um  of  (12).  Tha  phase  factor  ^  was  not  praaant  in  (25). 

An  aqulvalant  phaaa  factor  nay  ba  derived  directly  from  an  SPA  of  (24)  by 
assuming  that  u  is  large.  This  derivation  is  outlined  in  Appendix  B. 

A 

In  (29) .  B(b  )  is  given  by  (16)  with  6*6..  Tha  SPA  result  given 
n  n 

above  requires  |R(8q) |  to  ba  not  too  small  and  hence  that  6q  is  not  close 

to  Brave tar's  angle12  6^  (known  also  as  tha  angle  of  intromission) .  This 

angle  ia  defined  as  tha  solution  of  B(8g)  “  0;  inspection  of  (15).  gives 

ain26B  -  (M2  -  T2)/(l  -  T2).  If  T  <  M  <  1  or  if  T  >  M  >  1,  then  0B  exists. 

When  M  <  1.  all  a  ere  positive.  If  a  <0  (which  occurs  for 
n  A 

certain  M  >  1),  tha  signs  of  tha  ir/4  phases  in  (25)  and  (29)  should 
22 

be  reversed. 

C.  The  Axial  Bays 

Since  the  axial  rays  are  unfocused,  their  far-field  amplitudes  end  phases 

1  13 

can  be  computed  directly  from  ray  optics.  *  Figure  4  defines  the  distances 
needed  to  compute  the  phase  difference  Cq  between  the  n  ■  0  reflected  ray 
and  the  propagation  delay  from  C'  to  Q.  Using  our  (...)  notation  for  dis¬ 
tance  gives  Cq  »  k[(234)  -  (C'l)];  this  can  be  written 

-  2ka[l  -  coe(y/2)]  -  ka(l  -  cosy)  (32) 

by  using  the  geometric  results  8  ■  y/2  and  (C'l)  “  (56).  The  far-field 
pressure  due  to  each  axial  ray  is  p.[(a/2r)  exp(lkr)]f  where  the  reflected 

l  A 

ray  has  the  following  fora  function  when  ka  is  not  small: 

fQ  -  -R(0  -  y/2)eiCo  .  (33) 

The  minus  sign  results  from  our  definition  that  Bq.  (15)  describes  internal 

reflections.  Whan  y  is  small,  the  reflected  wave  appears  to  cone  from  AQ 

23 

which  (as  y  0)  is  a  distance  e/2  from  C. 
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The  form  functions  f  (n  ■  2.4.  .  .  .)  of  th«  other  axial  rays  mar  be 
— — ■  ■—  n 

found  by  computing  divergence  factors1  end  the  phase  shifts  due  to  peth  lengths 
end  foci.1  It  is  not  generally'  feasible  to  express  the  ray's  angle  of  Inci¬ 
dence  0  as  a  function  of  Y*  this  necessitates  either  the  numerical  solution 
of  transcendental  equations  or  the  use  of  approximations.  An  approximate 

description  of  f 2 (Y)  is  given  in  Appendix  C.  A  description  with  y  m  0  of 

14 

electromagnetic  axial  rays  from  a  dielectric  sphere  can  be  modified  to  give 
the  following  acoustic  result  when  ka  »  1 

|fft(Y-0)|  *  |K(n  -  K)"1  -  R(0)21|  .  (34) 

For  the  examples  to  be  described  In  this  paper  fg  and  f2  are  similar  in 
magnitude;  for  n  2  4.  |«.|  «  I  fjj|  so  that  approximations  for  those  fQ(Y) 
will  not  be  given  here. 

D.  The  Combined  Scattering  Amplitude 

The  above  results  may  be  combined  to  give  the  following  approximation 
for  the  pressure  amplitude  la  die  far  field 

p(r.Y)  ■  PjKe/lr)  exp(lkr)]  f,  r  »  ka2  (35) 

V 

f(Y)  *  f0  +  f«  I  g.  M  <  1  (36) 

u  n-3 

where  N  should  be  sufficiently  large  to  approximate  an  infinite  series. 

Th.  wnaUtaUan  tor  th.  tor.  faction  t  In  b«o  cho.«  .uch  thot  1 
for  geometric  reflection  from  a  fixed-rigid  sphere  of  the  same  else.  The 
salient  feature  of  (29)  is  that  ^(Y  "  0)|“  (ka)1^2  while  the  |fa|  do 
not  depend  on  ka.  Consequently  from  large  fluid  spheres 

will  be  dominated  bv  contribution  of  the  diffracted  glory  waves  provided  the 
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attenuation  of  sound  is  negligible.  Aa  ka  tha  glory  contributions  to 

(36)  diverge;  this  divergence,  a  consequence  of  axial  focusing,  was  also 

evident  in  the  purely  geometric  scattering;  e.g. ,  Sq.  (6). 

When  ka  is  not  large,  (36)  oust  be  modified  to  include  terms  due  to 

circumferential  waves.  Our  tests  of  (36)  had  ka  >  100  and,  as  a  consequence 

15  19  24 

of  the  largeness  of  ka,  circumferential  waves  *  and  associated  reaonances 

should  experience  significant  radiation  damning.  The  backs cattering  from 

spheres  due  to  circumferential  waves  will,  nevertheless,  be  assisted  by 

diffraction-limited  axial  focusing  as  is  the  case  for  electromagnetic 
1  3  12 

scattering.  *  *  Circumferential  wave  contributions  to  the  scattering  from 

25  26 

elastic  objects  are  significant  for  ka  as  large  as  *  200;  however,  their 

significance  to  the  scattering  from  large  fluid  objects  with  *  qq  is  not 
well  explored. 

10  13  24 

The  conventional  description  of  the  scattering  *  *  references  the 

phase  of  the  incident  wave  to  the  sphere's  center  C.  It  also  uses  distances 
and  angles  with  respect  to  C  so  that  r  becomes  (0Q)  and  Y  becomes  the 
polar  angle  of  Q  with  C  as  the  origin.  For  the  far- field  scattering  the 
form  of  (35)  is  retained  with  f  replaced  by  fg(Y)  ■  f(Y)  expUCg]  where 

■  -  ka(l  +  cosy)  (37) 

V 

is  the  negative  of  the  phase  shift  for  the  distance  (C'CS)  in  Fig.  4.  The 

Y-dependant  phase  shift  for  the  reflected  ray  becomes  ♦  Cc  •  -2ka  cos(y/2) 

which  agrees  with  the  conventional13  result;  for  the  glory  terms  it  becomes 

+  C  -  -k  (2a  +  (a  -  a) (1  -  cosy)].  (When  a  •  a,  tha  latter  result 
n  c  u  n 

follows  from  elementary  considerations.)  Since  the  modulus  of  the  form  func¬ 
tion  is  not  altered  by  the  transformation  from  C'  to  C,  it  is  possible  to 
compare  ) f |  from  (36)  directly  with  |fj  where  f#  is  the  exact  result  of 
the  partial-wave  theory.10'1®’3* 
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IV.  Discussion  of  Modal  sad  Comparison  with  Exact  Scattaring  from 
Spharas  with  M  <  1 


In  this  sactlon  wa  compara  Eq.  (36)  with  the  modulus  of  the  asset  form 
function  fg.  Tha  partial-wave  series  for  was  sumaad  by  using  the  com¬ 
puter  algorithm  described  in  tha  appendix  of  Be£.  18.  Tha  number  of  partial- 
waves  included  exceeded  ka  +  4.05  (ka)1^  to  ensure  adequate  convergence.^* 
This  algorithm  was  limited  to  the  case  of  equal  inner  and  outer  fluid  densities 
(p^  -  pg)  and  the  |f#|  presented  hare  and  In  Sec.  V  are  limited  to  this  ease. 

Tables  I  and  II  are  representative  of  model  results  for  form  function 
moduli  of  axial  and  glory  terms  and  for  tha  focal  circle  parameters.  These 
tables  should  be  examined  in  conjunction  with  Fig.  5-7.  In  each  of  the 
figures,  the  largest  H  was  chosen  to  be  a  power  of  2  such  that  an  additional 
doubling  of  the  largest  n  in  (36)  resulted  in  a  new  curve  (not  shown) 
which  differed  imperceptibly  from  the  dotted  curve. 

The  value  of  M  in  Fig.  5  was  selected  such  that  the  modulus  of  f^, 
the  strongest  of  the  omitted  axial  ray  amplitudes,  is  especially  smell  in  com¬ 
parison  with  |g3(Y  -  0)|.  The  — <n  of  Fig.  5a  is  that  (36),  with 

R  i  16,  gives  an  amplitude  which  is  nearly  identical  to  the  exact  result  for 
the  range  of  y  plotted.  This  confirms,  with  very  large  ka,  our  model  result 


for  the  amplitude  end  relative  phases  of  the  fQ  and  g^.  Vote  also  that  |f| 
for  y  >  0.2*  is  dominated  by  the  interference  of  gj,  fQ,  and  t^i  for 
y  <  0.2*,  however,  glory  terns  with  n  >  3  are  significant.  Figure  5b  con¬ 


firms  that  the  principal  features  of  | f  —  |  are  described  by  the  model.  Mote 

that  the  dip  la  |f#|  near  12.5*  is  largely  due  to  the  destructive  interference 

4*6  8 

of  axial  rays.  Though  we  have  previously  modeled  individual  contributions,  * 


Fig.  5  is  the  first  direct  confirmation  that 


logy  waves  can 
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Figures  6  and  7  confirm  again  chat  tha  principal  features  of  |f#| 
are  described  by  the  model.  In  Fig.  6a  and  7,  discrepancies  between  the  wo deled 
and  exact  scattering  are  evident,  especially  at  y  ■  0.  The  causes  of  these 
discrepancies  are  not  known.  Equation  (36)  has  also  been  tested  by  Including  a 
sua  of  glory  terms  with  n'  *  1  (which  have  n  >  5)  but  the  resulting  shift 
in  I f |  is  ouch  too  saall  to  account  for  the  discrepancies.  Though  the 
Individual  | fQ|  with  n  >.  4  are  saall,  it  la  plausible  that  the  coherent 
sua  of  omitted  axial  terns  could  account  for  significant  part  of  the  discre¬ 
pancy.  It  is  apparent  that  the  oaisslon  of  circuaferentlal  waves  froa  the  aodel 
is  acceptable  In  Fig.  5  and  6b;  this  oaisslon  could  account  for  soae  of  the 
discrepancies  evident  in  Fig.  7. 

Soae  noteworthy  features  of  the  modeled  scattering  are:  (1)  the  width 

of  the  backward  peak  of  the  scattering  is  roughly  *  l/(ka)  but  the  details 

of  the  structure  are  highly  dependent  on  ka  according  the  interference  of  the 

terns  in  (36);  (11)  there  is  a  tendency  for  the  width  of  the  peak  to  increase  with 

1/2 

decreasing  b^  and  hence,  decreasing  M;  (ill)  though  the  Ig^Y  ”  0) |  «  (ka)  , 

| f |  is  not  <*  (ka)  '  due  to  the  ka  dependence  of  the  interference  between 
the  g^  and  the  ka- in dependent  f^;  (lv)  fora  function  moduli  can  exceed 
(e.g.  Fig.  6a)  or  be  close  to  (e.g.  Fig.  6b,  7b)  unity  which  is  the  geometric 
result  for  reflection  froa  a  fixed-rigid  sphere  of  the  saae  else.  This  enhance¬ 
ment  is  a  Manifestation  .of  diffraction-limited  axial  focusing  of  the  back- 
scattering.  It  is  possible,  however,  for  the  various  glory  tens  to  Interfere 
eo  as  to  produce  a  ainiss»  in  |f  |  for  y  »  0;  this  is  evident  la  plots  of 
| f  |  and  |fg|  for  M  “  0.5  with  ka  »  1000  which  are  not  shewn  here. 

Soae  aspects  of  the  convergence  of  the  series  la  (36)  sould  be  noted. 

As  a  ♦  •,  ba  ♦  M  and  no  ■*  where 


n.  -  2k*  [(1  -  co*0c)  +  M(ir  -  8C)]  .  (38) 

This  1*  th«  propagation  phase  shift  a  ray  would  hav*  if  It  an tar ad  and  exited 
froa  tha  sphara  at  tha  critical  angla  6£  and  It  travalad  a  circumferential 
path  with  a  phaaa  velocity  of  c^.  Numerical  taata  suggest  that  It  Is  suf¬ 
ficient  to  ternlnate  the  series  wheri”  (n  -  n  )  <  1  radian.  The  omitted  tens 

•  n 

tend  to  cancel  because  of  the  alternating  sign  of  B(bQ)  and  the  periodicity 

of  the  phase  factor  exp(iy  )  associated  with  the  crossing  of  Internal  foci. 

n 

As  n  Increases,  |gQ |  decreases  due  to  decreasing  |B(bQ)|  and  Increasing 
q^;  see  Tables  I  end  II.  When  i  pQ,  the  B(ba)  are  changed  but  the  geo¬ 
metric  parameters  are  not. 

It  nay  be  possible  to  arrive  at  (29)  for  gQ  froa  the  asymptotic 
evaluation  (at  large  ka)  of  the  exact  partial-wave  sua.  Glory  ray  amplitude* 
should  be  dascrlbable  by  saddle-point  contribution  to  a  contour  integral  froa 

3 

the  Watson  transformation  of  the  exact  sua.  Our  us*  of  the  r’ysical-optlcs 

approximation  in  Sec.  II  and  III,  though  less  direct,  manifests  the  physical 

significance  of  the  paraaeters  a  ,  b  ,  q_  ,  <p  ,  and  u.  Furthermore,  our 

H  a  a  u  a 

approach  bypasses  the  difficulties  with  the  asymptotic  method  noted  In  Ref.  13 

and  It  may  be  extended  to  the  near  field  by  replacing  (19)  by  a  Fresnel 
10*21 

transform.  It  aey  be  important  that  p±  and  PQ  be  similar  In  aagnl- 

-  24 

tud*  for  otherwise  resonances  (as  In  s  gas  bubble)  aey  be  significant.  The 
lower  limit  on  ka  for  which  (36)  Is  applicable  Is  not  known;  the  physical- 
optics  aathod  was  found  to  be  useful,  depending  on  M,  for  the  description  of 
near  critical-angle  scattering  when  ka  i  25. 

The  physical-optics  model  aay  be  used  to  approximate  the  scattering  of 
ton*  bursts  where  k  Is  obtained  from  the  average  frequency  of  the  incoming 
burst.  The  time  delays  of  dlscrets  schoes  follows  from  the  propagation  phase 
shift*  n  ,  9  ,  end  c  .  Shapes  of  discrete  glory  echoes  will  differ  froa 

II  Tn  Q 
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eh*  Incoming  bursts^  due  to  th«  k- independent  phase  shifts*  tbs  yQ  and 
ir/4  tans  in  (29) ;  tha  scattarad  burst  is  ralatad  to  tha  incoming  signal 
through  superposition  of  tias-shif tad  incoming  and  Hilbart-transformad 
signals. ^ 

An  extensive  test  of  (36)  could  be  carried  out  bp  repeating  the  com¬ 
parisons  in  Figs.  5*7  with  several  smaller  ka.  This  would  be  inefficient  to 
do  at  present  since  the  |f#|  and  |f|  curves  are  plotted  using  separate 
computer  systems.  Figures  5-7  were  made  by  overlaying  and  tracing  the  curves. 

V.  Spheres  with  M  >  1  and  Combined  Kalnbov  and  Glory  Scattering 

In  this  section  we  summarise  the  results  of  exact  calculations  of 

backscattering  from  large  fluid  spheres  with  certain  M  >  1  and  we  model  the 

enhancement  of  backscattering  due  to  rainbow  rays.  As  described  in  Sec.  I, 

the  class  of  glory  rays  with  n*  »  0  and  n  >.  3  are  not  limited  to  M  <,  1; 

they  also  exist  for  1  <  M  <  M1  where  the  upper  bounds  are  given  by  (4)  and 

n 

(5)  for  n  ■  3  and  4.  Fora  functions  for  this  class  of  rays  ere  given  by  (29) 

except  when  the  ray  has  n  internal  chords  and  M  3  M7 .  There  is  an 

additional  caveat  to  be  noted: *  whan  If '  <  M  <  M'  there  is  second  class  of 

n  n 

rays  having  a  internal  chords  which  also  cross  the  axis  once.  Here 

M"  "  csc[(n-l)ir/2n]  is  the  value  of  M  for  which  (1)  and  (3)  give  0  -  6  •  ir/2. 

D 

When  n  -  3,  the  angle  of  incidence  of  this  second  backs cattered  ray  is  given 

by  (17)  with  (I*3  +  ir)/3  replaced  by  ir  +  (r3/3). 

The  following  coaparlson  of  ray  properties  for  the  two  classes  (each 

with  n  chords)*  will  facilitate  a  description  of  the  unusual  backscattering 

properties  of  spheres  having  M  3  M'.  Let  b  ■  sin8  and  a. (from  Eq.  (10)] 

n  s  n  a 

denote  the  focal  parameters  for  the  new  class  of  ray  while  b  and  a.  denote 

is  n 

those  for  the  original  type.  Vlth  M*J  <  M  <  m£,  the  parameters  obey  the 
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following  inequalities :  b  <  b  <  a,  o  >  a,  and  (for  M  not  too  close  to 

n  n  n 

M* )  a  >0.  As  M  M  ,  Eqs.  (1),  (3),  (10),  and  (17)  lead  to  the  results: 
u  n  n 

b„  ■+  b  ,  nT  ■*  1,  o  ♦  +»,  5  -*■  -»  and  q_  -*■  1.  The  divergence  of  the  die- 

tances  to  the  focal  circles  gives  rise  to  the  erroneous  prediction  by  (31) 

that  EqI  and  hence  gQ  also  diverge  for  each  of  these  rays.  The  present 

physical-optics  approximation  fails  because  the  glory  wave  is  no  longer  toroidal 

as  required  by  (13).  Numerical  computations  applied  to  the  n  ■  3  rey  and  the 

case  M  •  give  the  following  limit  as  s  ■+•  b^:  tn(s)  -  (•  ~  b^)5  -*■ 

_2 

Ajk a  where  the  dimenaionaless  constant  A^  *  -26.6.  Hence  the  wavefront  is 

cubic  which  is  characteristic  of  a  "rainbow"  or  "stationary"  ray.  Figure  8 
shows  this  ray. 

1  7  13  23 

It  is  well  known  that  scattering  is  enhanced  *  ’  ’  In  the  vicinity 
of  a  rainbow  ray.  Rainbow  rays  have  dy/dS  •  1  -  (dB/d6)  ”  0.  This  is  equiva- 

C  1  ^ 

lent  to  the  condition  |dbj/d$|  •*  ®  in  (6).  These  conditions  require  that  * 

sin20  -  (n2  -  M2)/(p2  -  1).  We  have  verified5  that  the  glory  rays  with 

M  «  and  n  -  3  and  4  satisfy  this  requirement;  geometrical  considerations 

suggest  that  it  is  also  mat  when  n  >  4.  Consequently  when  M  *  M' ,  the  back- 

n 

scattering  is  doubly  enhanced  relative  to  that  due  to  axial  rays:  once  due  to 

axial  focusing  and  once  due  to  the  rainbow  caustic. 

He  have  developed  a  physical-optics  approximation  for  gQ  for  the 

special  case  M  »  M' .  The  principal  change  of  the  approxlmetion  described  in 
n 

Sec.  I1IA  is  to  replace  (24)  by 

I  -  £.  J0(.).lfl«k(-b")5'*J  d.  .  bnJ0(»n)It  (3») 

where  the  SPA  has  been  used  to  remove  from  the  integral,  that  part  of  the 
integrand  which  is  assumed  to  be  slowly  varying  near  s  -  bQ.  The  remaining 
integral  is 
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r  #iAnk(8-bn)3/a2 

Jo 


(40) 


la  the  analysis  which  follows,  we  assume  chat  Aq  <  0,  which  is  Che  case  for 
a  •  3.  It  can  be  shown  chac  Che  general  form  of  Che  final  result,  Eq.  (42), 

does  noc  depend  on  the  sign  of  Aq.  Changing  Che  integration  variable  to 

-2  1/3 

4  -  (s  -  bQ)  (~3AQa  )  reduces  1^  Co  an  incomplete  Airy  integral  in 

-2  1/3 

which  the  lower  limit  of  integration  is  A  ■  -b  (-3A  a  )  .  As 

a  n  n 

1/3 

k  4q  ■*  -0°,  the  incomplete  integral  has  an  asymptotic  approximation  [Eq.  (10) 
of  Ref.  28]  in  terms  of  the  complete  Airy  function  Ai  which  gives 


2irAi(0) 


+  i<r> 

n 


2  e 


-i(A'/3) (b-/ a) 


^2/3 


(41) 


where  Ai(0)  =  0.35503  and  A*  -  -3kaA  .  These  approximations  give  the  fol- 

n 

lowing  normalized  form  function  in  place  of  gQ  when  M  “ 


hn  *  -‘(k*)273  E«  ■Oa)J0<V1In°  +  “”I 


(42) 


E  -  Air  Ai(0)  b  /[a(-3A)1/3] 
/lx i  a  a 


(43) 


where  we  have  used  the  geometric  result  chat  qQ  •»  1  and  the  second  term  in 

(41)  has  been  omitted  because  of  its  small  magnitude  in  cases  to  be  considered 

2/3 

here.  It  is  to  be  expected  that  g.  *  (ka)  because  the  enhancement  of  the 

'in 

1  3  23  1/6 

diffracted  amplitude  of  a  nonbackscatterad  rainbow  *  *  is  «  (ka) 

The  results  of  this  model  when  n  ■  3  era  compared  with  |fe|  in 
Fig.  9.  Away  from  y  •  0,  a  complete  approximation  for  g^  may  contain  a 
y-dapendant  phase  shift  similar  to  the  factor  in  gQ.  Since  this  shift 

has  not  been  determined,  our  comparison  is  limited  to  comparing  ) g^ I  with 
|fcl  in  cases  where  | f# |  »  | f 2  j  (which  is  now  the  leading  axial  ray  amplitude. 


22 

sac  Table  III).  Figure  9  demonstrates  that  (42)  describes  the  main  features  of 

the  scattering  when  ka  ■  1000  but  that  (42)  is  relatively  incomplete  when 

ka  ■  100.  This  conclusion  is  also  supported  from  our  computations  of 

jf^Cy  “  0) |  for  several  ka  within  ±10  of  1000  and  100.  In  each  case  |f#| 

varied  aperlodically  with  ka.  Hear  ka  “  1000  the  extrema  were  typically 

between  17.4  and  20.2;  near  ka  *  100,  they  were  typically  between  3.8  and  6.8. 

These  are  to  be  compared  with  of  18.9  and  4.1  for  ka  of  1000  and  100, 

respectively.  That  there  should  be  r.gnif leant  corrections  to  the  physical- 

optics  approximation  when  -ka  -  100  is  to  be  expected  from  proximity  of 

-  55.7°  to  the  internal  critical  angle  ■  arcsinOl^-1)  -  57.9*.  The 

use  of  plane-surface  reflection  coefficients  to  compute  B(b,)  fails  as  v  -*■  v 

J  c 

because  tunneling,  surface  waves,  and  resonances  make  curvature  an  important 
1  3 

consideration.  ’ 

Table  III  includes  exact  results  with  other  values  of  H.  For  each 
case  | |  was  computed  for  a  range  of  y  sufficient  to  include  several  of 
the  backward  diffraction  maxima;  its  peak  occurred  at  y  •  0  in  each  case. 

The  relative  magnitudes  are  consistent  with  the  considerations  given  at  the 
beginning  of  this  section.  The  rainbow  condition  is  not  met  when  M  ■  1.1  and 
the  backscattering  is  weaker  than  for  M^.  The  backscatterlng  is  doubly 
enhanced  when  M  ■  Mj.  It  is  small  when  M  ■  1.25  with  ka  -  1000  because 
there  is  no  glory  ray  having  only  a  few  chords  (irrespective  of  o')  and  circum¬ 
ferential  waves  experience  significant  radiation  damping. 

The  combined  rainbow-glory  enhancement  of  backscattering  is  not  limited 

3 

to  the  precise  condition  M  ■  M^.  Due  to  diffraction,  a  rainbow  will  influence 
the  backscattering  when  dy/d8  ■  0  in  the  vicinity  of  y  •  0.  Consequently, 
beckscattering  will  be  enhanced  for  M  in  some  range  near  which  narrows 
for  increasing  ka.  There  are  other  cases  of  glory  enhanced  backscattering 
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outside  of  the  range  0  <  M  -  considered  in  this  paper.  The  most  useful  of 

these  may  be  that  due  to  the  n  »  2  glory  ray  for  Jl  <  M  <  2. 

29 

Our  Mie  theoretic  computations  of  the  backscatterlng  of  light  from 
dielectric  spheres  also  demonstrate  a  combined  rainbow-glory  enhancement  for 


optical  refractive  Indices  of 


and  M^. 


VI.  Applications  to  Underwater  Acoustics 

The  models  developed  in  Sec.  Ill  and  V  for  glory  and  rainbow-enhanced 
glory  backscatterlng  are  more  general  than  the  case  of  a  purely  fluid  sphere 
provided  focal-circle  parameters,  attenuation,  and  coefficients  of  reflection 
and  transmission  are  properly  modeled.  For  example,  Eq.  (29)  removes  the 
Y  -*■  0  divergence  of  amplitudes  backscattered  from  solid  spheres  present  in  the 

geometric  model  of  Ref.  9;  subtleties  of  this  application  will  be  described 

XX  30*32 

elsewhere.  Spheres  have  been  used  as  calibration  targets  for  sonar  systems'3 

because  of  their  symmetry;  their  response  is  asserted  to  be  more  uniform  (in 

regard  to  variations  in  the  direction  of  the  incoming  wave)  than  the  triplane 
31  32 

reflector.  ’  In  this  section  we  comment  on  the  design  of  practical  spheres 

which  could  be  made  to  exhibit  glory  enhanced  backscatterlng. 

Liquid  spheres  with  ka  -  100  are  too  large  for  surface  tension  to 

ensure  sphericity.  For  example,  with  c0  ■  1.53  m/s  (sea  water)  and 

u>/2ir  ■  100  kHz,  the  radius  =  24  cm  when  ka  ■  100.  The  target  liquid  may  be 

31  32 

contained  in  a  thin  elastic  shell.  ’  For  ease  of  transport,  the  shell's 

31 

interior  may  be  left  unfilled  until  submersion.  In  the  discussion  which 
follows,  shear  and  longitudinal  sound  velocities  of  the  bulk  shell  material 
will  be  denoted  by  c^  and  c^  and  the  shell's  thickness  will  be  denoted  by 
h;  p^,  a,  and  c^  are  the  density,  radius,  and  sound  velocities  of  the  inner 


liquid  and  M  ”  cQ/ci 
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For  a  large  sphere,  glory  rays  illustrated  in  Fig.  1  and  8  exist  In 

the  h  «  a  Halt.  The  main  difficulty  in  using  Eq.  (29)  and  (42)  to  estimate 

the  scattering  Is  in  the  evaluation  of  B(bQ)  for  the  tranaalssion  through 

and  reflection  from  the  shell.  Proper  choice  of  h  could  simplify  the 

analysis  by  facilitating  the  use  of  the  thin-plate  approximation  for  the 

17  19 

transmission  and  reflection  coefficienta  of  flat  plates.  *  For  larger  h, 

the  analysis  is  complicated  since  both  c.  and  c  typically  exceed  c  and 

m  8  O 

24 

c^.  Furthermore,  Lamb  modes  are  launched  in  the  shell  at  certain  angles  of 
26 

incidence.  A  plausible  choice  to  facilitate  transmission  into  the  sphere's 
interior  is  to  select  c#  somewhat  smaller  than  cq  and  c^.  (For  example, 
Lexan  is  a  atrong  material  having  c^  -  910  m/s.)  In  the  case  of  rainbow- 
enhanced  glory.  Fig.  8,  the  calculation  of  B(b3)  may  be  especially  complicated 
due  to  the  largeness  of  8^  and  v^. 

A  computational  demonstration  of  glory-ray  enhanced  backseat taring  may 
be  best  achieved  by  computing  the  exact  scattering  from  thin,  liquid-filled 
apherlcal  shells  by  extending  the  computations  in  Ref.  33  to  ka  Z  100  and 
certain  MM.  The  calculations  should  include  a  range  of  Y«  The  exact 

2/3 

calculation  would  be  most  interesting  in  the  case  M  ■  due  to  the  (ka) 

enhancement  factor  in  (42).  The  possibility  of  enhanced  backscattaring  from 

liquid-filled  shells  with  this  M  has  been  previously  overlooked.  (For 

example,  if  it  is  present,  the  interpolated  target  strength  versus  M  curve 

in  Fig.  3  of  Ref.  31  contains  serious  errors.)  Circumferential  waves,  which 

26 

are  known  to  influence  backscattaring  from  large  cyllnderical  shells,  will 
give  rise  to  axially  focused  backscattaring  from  spherical  shells.  The  relative 
importance  of  ray-optical  and  circumferential  returns  may  be  evaluated  by  dis¬ 
playing  the  scattering  of  a  short  tone  burst  as  in  Ref.  33.  Certain  inhomo- 

34 

geneous  spheres  exhibit  glory  rays  and  should  have  enhanced  backscattaring. 
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Appendix  A:  Curvature  and  Spreading  of  Glory  Wavefronts 

The  wavefront  is  described  by  a  curve  of  constant  propagation  phase 
delay  n>  Its  curvature  at  tha  exit  plane  is  1/a  “  k  A(d  n/dsA)/[l  +  (k  dn/ds)  ]  ' 
where  s  is  the  distance  free  C';  n  and  s  are  given  by  (7)  through  (9)  as 
functions  of  8,  8.  and  v,  but  could*  in  principle,  be  written  as  functions  of 
6  only,  by  making  use  of  (1)  and  (3).  Me  then  have 

1  .  (da/d0)(d2n/d02)  -  (dn/d6)  (d2s/d02) 

°  k[(ds/d8)2  ♦  (k*1  dn/dS)2]3/2 

where  the  differential  operator  is 


(Al) 


jL  m  _5_  4.  (iLi)  JL  4.  &L  4.  _L 

d8  30  'd0'  3v  l3v  d0  30;  36 


-  3/36  +  T  3/3v  +  (2nx  -  1)  3/36 


(A2) 


(A3) 


where  (A3)  uses  (1)  and  (3)  and  tha  definition  T  ■  dv/d8  »  tanv/tan6.  The 
first  derivatives  in  (Al)  are 

(ka)”1  -2  (l-nT)(l-cos6)sec(0-8)tan(0-6)  +  (l-2nT)(sin0-sin6sec(0-6))  ,  (A4) 

a”1-  -  2(nT-l)(l-coa6)eac2(0-8)  +  C  nT-l)(cos8-sin6tan(8-8))  .  (A5) 


Tha  second  derivatives  are  somewhat  longer,  but  straightforward  to  calculate. 
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The  Interesting  cases  are  the  glory  rays ,  occurring  when  6*6;  this  condition 

will  be  denoted  by  e  subscript  n  for  the  glory  ray  with  n  chords.  Since 

(dn/dd)  -  0,  (Al)  becomes  k/a  -  [(d2n/d26)/(ds/d6)2l  .  The  derivatives  at 
non 

the  glory  condition  are  (ds/d6)  ■  a[2(nf-l)  +  cos6_]  and  (d2n/d02)_  • 

n  n  n 

2k(nT-l)(ds/d6)  .  Hence  the  radius  of  curvature  of  the  wavefront  is  given  by 
n 

(10b) . 

From  (9)  and  (11),  the  spread  of  the  wavefront  aey  be  written  as 


.  11* 

0,6*6.. 


a in 8  -  [sing  -  (1  -  cosS)tan(6  -  B) ] 

a 

sin6  -  sin€ 
n 


(d/d9)[sin0n  -  sinB  +  (1  -  cosB)tan(6  -  6)1 

0,6*6..  (d/d8) [sin0  -  sin0] 

n  n 


where  (A7)  follows  from  L' Hospital 'a  rule.  Application  of  (A3)  gives  an 
expression  which  reduces  to  (12). 

We  have  verified  (10b),  for  several  values  of  M  +  H’ ,  with  direct 

II 

2 

numerical  computations  demonstrating  that  as  s  ■*  b  ,  [n(s)  -  n  ]/(s  -  b  )  •» 

n  u  q 

k(2a  )  Furthermore,  determinations  of  focal  circle  locations  by  direct  ray 

U 

tracing  (from  large  versions  of  Fig.  1)  are  in  agreement  with  (10b).  Our 
result  for  has  also  been  verified  by  numerical  evaluation  of  the  limit 
in  (11).  These  tests  were  merited  because  quantities  equivalent  to  a  end 

a 

qQ  are  given  by  incorrect  expressions  in  Appendix  II  of  Kef.  23.  Those 
expressions  were  erroneous  due  to  incorrect  formulation  of  total  derivatives; 
they  happen  to  give  the  correct  ratio  1^1/^* 

Appendix  B:  Angle-Dependent  Phase  Shift  Via  the  Method  of  Stationary  Phase 
The  purpose  of  this  appendix  is  to  demonstrate  that  a  modified  SPA  of 
(24)  yields  a  phase  shift  equivalent  to  (30)  and  to  give  insight  into  the  cause 


27 


o£  that  shift.  The  derivation  which  follows  is  limited  to  cases  where 
uQ  »  1.  In  (24),  express  JQ(u)  using  Henkel  functions  of  the  first  and 
second  kinds.22  Jq(u)  »  (u)  +  Hg2^(u))  and  define  Zj(u)  •  Hg^  (u)expQlu) 

with  j  -  1  end  2  (here  and  below)  for  the  upper  and  lower  sign,  respectively. 
Then  (24)  is  a  sun  of  integrals  of  the  fora 

Ij  •  %f.  Zj  exp[ikJi(s-bQ)2  a*  ±  iuj  ds  (Bl) 

where  the  stationary  phase  points  of  the  coup lex  exponentials  are  at  *  "  *j 

with  s.  ■  b  r  et  T,  Froa  the  asyaptotlc  forms22  of  the  H^(u)  as  u  ■*>  «®, 
j  n  T  n  u 

it  is  evident  that  sZ^  is  a  slowly  varying  function  of  s  near  s  •  s^  pro¬ 
vided  u  »  1  and  jo  |T  «  b  .  Consequently,  in  the  SPA  of  (Bl),  it  is 
appropriate  to  treat  the  coaplex  exponential,  which  differs  froa  that  of  (24), 
as  the  function  which  oscillates  rapidly  when  s  +  s^.  Approximating  1^  and 
I2  by  this  procedure  gives  the  following  sum: 

I  *  bn(2x|on|/k)1/2  J0(un)e1(*“  +  ff/4)  (B2) 

where  terms  of  order  aT/b  relative  to  unity  have  been  neglected  and 

n  n 

<r '  -  -  >ska J2.  Replacing  (25)  by  (B2)  leads  directly  to  (29)  with  <p  replaced 
*  u  n  1  q 

by  (p\  By  Inspection,  when  T  «  1,  <pn  and  are  identical  up  to  terms 

4 

«r  .  They  give  nearly  identical  results  for  the  conditions  under  which  (29) 
was  tested  in  Sec.  IV.  The  derivation  of  ©',  however,  assumes  that  both 

1  a 

u  »  1  and  |a  | r  «  b  and  hence  that  (kb  )*2  «  T  «  b  /|a„|. 

The  shift  of  the  stationary  phase  points  by  t  aQr  has  the  following 
physical  interpretation.  Consider  the  locations  of  the  effective  areas  (or 
Fresnel  sones2)  of  the  toroidal  wavefront  which  contribute  to  the  scattering  to 
Q.  Whan  r  #  0,  these  tones  are  cantered  on  points  with  tp  ■  £  snd  ip  ■  ir  +  £ 
where  £  (Fig.  3)  is  the  aslnuth  of  Q.  The  centers  of  these  cones  are  shifted 
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away  from  b  by  i  a  T.  This  shift  of  tha  affaetiva  araaa  of  eha  toroidal 

u  n 

wava  laads  to  tha  phaaa  shift  <p\ 

Appandix  C:  Aaplltuda  of  tha  One-Bounca  Axial  Bay 

Tha  purposa  of  this  appandix  Is  to  dascrlba  an  approximation  for 
f2(y).  Tha  geometry  of  tha  assoelatad  ray  Is  shown  In  Tig.  C-l.  Lat  C2 

danota  tha  phaaa  dlffaranca  batwaan  tha  n  •  2  ray  and  tha  propagation  dalay 

from  C'  to  Q.  Inspactlon  of  Fig.  C-l  givaa  C2  “  k[(12)  +  M(234)  +  (45)  -  (C'6)]  • 
2ka[(l  -  cos6)  +  2M  cosv]  -  ka(l  -  cosy)*  We  axprasa  (2  ai  a  function  of  y  by 
using  an  approximation  for  €.  As  y  -►  0,  Eq.  (1)  bacoaas  8  *  Mv  provldad  |M  -  l| 

Is  not  larga.  Eliminating  v  from  tha  axact  axprasslon  y  ■  2(2v  -  8)  givaa 

8  »  My/(4  -  2M)  (Cl) 

From  this  astlaata  of  8,  tha  complata  Eq.  (1)  Is  usad  to  obtain  v  and  c2 
Is  found  via  tha  axprasslon  givaa  abova.  This  approxlaatlon  for  C2  ***  com- 
parad  with  tha  axact  rasult  (which  nay  ba  conputad  as  a  function  of  8).  Tha 

arror  Is  nagliglbla  for  tha  ranga  of  ka,  y,  and  M  of  tha  computations  in 

Sac.  IV.  Tha  approximation  for  f2  is 

f2  *  B(8)(l  -  R(0)2][M/(2  -  M)]ai(C2"W)  (C2) 

whara  tha  -it  phaaa  tarm  rasults  from  tha  crossing  of  two  foci.  Tha  factor 
M/(2  -  M)  Is  tha  dlvarganca  factor^’^^*^^  approprlata  for  y  •  0.  Tasts 
indlcata  that  our  approxlaatlon  of  this  factor  by  a  constant  lntroducas  a 
nagliglbla  arror  In  Fig.  5-7.  Tha  raflaction  coefficient  was  computed  via  (15) 
and  (Cl). 
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TABLE  II.  Focal  circle  parameters  for  M  -  0.94  (upper  group)  and 
0.6  (lower  group).  B^)  la  shown  for  che  case  Pi  -  pQ. 


n 

0Q(deg) 

V* 

E 

n 

B(bn) 

3 

40.24 

0.646 

1.162 

7.2 

1.30 

3.0  E-3 

4 

53.08 

0.800 

1.078 

13.9 

1.12 

-9.1  E-4 

8 

65.62 

0.911 

1.016 

65.1 

0.57 

-1.2  E-4 

16 

68.91 

0. 933 

1.004 

277.1 

0.28 

-4.0  E-5 

3 

21.20 

0.362 

1.096 

11.4 

0.56 

9.3  E-2 

4 

28.24 

0.473 

1.052 

20.4 

0.54 

-5.7  E-2 

8 

34.76 

0.570 

1.012 

84.2 

0.31 

-2.3  E-2 

16 

36.33 

0.593 

1.003 

340.3 

0.16 

-1.1  E-2 

TABLE  III 

selected 

Exact  and  axial 

M  >  1  when  p.  •  p 
1  0 

form  function  moduli  at  y  • 

• 

1 

0  for 

M 

|f€(ka  -  100)j 

|f#(ka  -  1000) | 

lfol 

1**1 

2.57 

7.12 

0.042 

0.051 

1.10 

1.97 

1.99 

0.048 

0.058 

*3 

6.50 

19.60 

0.083 

0.118 

1.25 


1.39 


0.26 


0.1U 


0.183 
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Figure  Captions 

Fig.  1.  Backscattered  rays  from  a  sphere  with  M  *  0.6.  The  center  of  the 
sphere  is  C  and  the  figure  may  be  rotated  about  the  CC'  axis. 

Fig.  2.  Path  of  a  ray  (dashed  line)  as  it  leaves  the  sphere. 

Fig.  3.  Angles  and  distances  needed  to  describe  a  point  (x* ,y')  in  the  exit 
plane  and  the  observation  point  Q.  The  z  axis  is  the  extension  (toward 
the  source)  of  the  CC'  axis. 

Fig.  4.  Distances  needed  to  describe  the  phase  of  the  reflected  wave. 

Fig.  5.  Comparison  of  exact  and  model  form  functions  of  spheres  with  M  ■  0.94. 
The  curves  labeled  N  give  |f|  from  Eq.  (36)  evaluated  with  the  indicated  N. 
In  (a),  the  model  result  (dotted  curve)  is  nearly  identical  to  the  exact  result 
(solid  curve). 

Fig.  6.  Comparison  as  in  Fig.  5  but  with  M  *  0.6. 

Fig.  7.  Comparisons  as  in  Fig.  5  but  for  (a)  M  »  0.8,  and  (b)  M  -  0.5,  with 
ka  “  100  in  both  cases. 

Fig.  8.  Combined  rainbow-glory  ray  for  a  sphere  with  M  «  of  Eq.  (4). 

The  distance  scales  in  the  sketch  of  the  cubic  wsvefront  have  been  enlarged  for 
enhanced  visibility. 

Fig.  9.  Comparison  of  exact  (solid  curves)  and  modeled  (dashed  curves)  form 
functions  of  spheres  with  ka  »  1000  (upper  group)  and  100  (lower  group).  The 
model  results  are  from  Eq.  (42)  which  describes  only  the  rainbow-enhanced  glory 
ray.  The  ray  is  the  n  •  3  rsy  in  Fig.  8  which  has  b^  =  0.9752  a. 


y  (degrees) 


CUBIC  WAVE  FRONT 
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Submitted  to  the  104th  Meeting  of  the  Acoustical  Society  of  America  ! 

9-12  November  1982  ] 

Observation  of  the  acoustic  glory:  Scattering  from  an  elastic  sphere  in  near 

backward  directions.  Philip  L.  Marston,,  Timothy  J.  B.  Hanson,a^and  Kevin  L.  j 

i 

Williams  (Oept.  of  Physics.  Washington  State  University,  Pullman,  WA  99164) 

We  have  measured  the  scattering  for  small  angles  y  (relative  to  the 
backward  axis)  from  a  fused  silica  sphere  of  radius  <*■  %  52  mm.  Tone  bursts 
In  water  corresponding  to  k a  %  450  were  incident  on  the  sphere;  their  short 
duration  permitted  glory  and  axial  returns  to  be  separated  In  time.  The  y 
for  the  probe  hydrophone  was  scanned  to  test  a  model  [P .  L -  Marston  and 
L.  Flax,  J.  Acoust.  Soc.  Am.  Suppl.  68,  S81  (1980)]  of  diffractive  effects  on 
backward  axial  focusing.  Observations  tend  to  support  the  model  as  adapted 
to  fused  silica:  (1)  from  the  arrival  time,  the  strongest  echo  is  evidently 
due  to  the  4-chord  shear  glory  ray;  (2)  Its  amplitude  is  «  JQ(kbsiny)  where 
b  Is  the  calculated  glory  circle  radius;  (3)  Its  amplitude  at  y  *  0,  though 
slightly  smaller  than  predicted,  exceeds  that  of  the  first  axial  reflection; 
and  (4)  the  times,  amplitudes,  and  y  dependences  of  other  echos  are  correlated 
to  predictions.  The  first  null  of  the  strongest  echo  occurs  at  y  =  1°. 

Consequently,  we  demonstrate  for  the  first  time  the  diffraction  limited  backward 
focusing  of  echos  from  a  sphere.  [Work  supported  by  ONR.  Marston  is  an 
Alfred  P.  Sloan  Research  Fellow] 
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